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Abstract 

The problems with the Zjy symmetry breaking in the induced QCD are analyzed. 
We compute the Wilson loops in the strong coupling phase, but we do not find the 
Zn symmetry breaking, for arbitrary potential. We suggest to bypass this problem 
by adding to the model a heavy fermion field in a fundamental representation of 
SU(N). Remarkably, the model still can be solved exactly by the Rieman-Hilbert 
method, for arbitrary number Nj of flavors. At Nj <C N — > oo there is a new 
regime, with two vacuum densities. The Zn symmetry breaking density satisfies 
the linear integral equation, with the kernel, depending upon the old density. The 
symmetry breaking requires certain eigenvalue condition, which takes some extra 
parameter adjustment of the scalar potential. 
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1 Introduction 



1.1 Critical phenomena in induced QCD 

The first few months of investigation of induced QCD [0, f| |, |, |, |, |, |, [TO], [TT|, 
revealed some nice features of this new lattice gauge model, as well as some problems. 
The functional integral of this model reads 

Z = H [D$ 9 exp-Ntx [U($ x )] J] I[&x,$ v ], (1.1) 

x J <xy> 

with the Itzykson-Zuber integral 
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where the scalar field <& x is the N ® N hermitean matrix and the gauge field Q xy is the 
N (g) N unitary matrix. For the purposes of induction of QCD at large scales there is 
nothing wrong with the U(l) subgroup of U(N), therefore we leave the trace of $ finite, 
and det Q ^ 1 to simplify the equations. 

The first nice feature is the unique possibility to solve the model analytically at iV = oo 
by reducing it to the nonlinear integral equation (master field equation) for the vacuum 
density p(X) of eigenvalues of the scalar field [@], and the linear one for the fluctuation 
5px(X) I- 

The dimension D of the lattice enters as a parameter in these equations. The critical 
phenomena are quite rich: there are infinite number of fixed points p a (X) = \X\ a , with 
different critical indices a = 2n + 1 ± ~ arccos 3£ ^ 2 . Explicit massive solutions, interpolat- 
ing between pairs of the fixed points p a (X) and p2k-a(X) were found in the last paper JTT| . 
In particular, the masses of scalar excitations scale as 

m phv .<xB;6=-+ ^ _ fc) , (1.3) 

where B is the coeficient at the perturbation 

p (X) = \X\ a + B\X\ 2k - a ; a>k. (1.4) 



At the critical point this coefficient vanishes, which is needed to suppress the oper- 
ator of the lower scaling dimension. One would expect this B to vanish as some lin- 
ear superposition of deviations of the parameters of the scalar field potential Z7(3>) = 
\rriQ^ 2 + Ao$ 4 + . . . from the critical values. As it was discussed in flT| , the multicritical 
point of the type a, k, m is realized at a > 2m + 1 when the terms O ($ 2m ) in the scalar 
potential become relevant in the wave equation. 

The model, therefore, has a nontrivial thermodynamics with calculable indexes in 
arbitrary dimension. The forbidden interval is \ < D < 1, and the solutions at D > 1 



1 



show no pathologies, contrary to the string models. The string models can be regarded [12 
as linear realizations of the gauge symmetry in the same lattice gauge model, with the 
unitary measure replaced by the Gaussian one, like in the Weingarten model. This leads 
to the singularity at D = 1 with the forbidden interval D > 1. 

The difference between the unitary and the Gaussian measure is known to produce 
an extra structure at a random surface in the strong coupling expansion. In the weak 
coupling expansion, the usual perturbative QCD is hopefully recovered in the vicinity of 
the unit element, which arises as the classical solution for the unitary measure, but has 
no special role in case of the Gaussian one. 

In a way, this field theory of the field p x {ty, with extra continuum coordinate A 
represents the theory of extended objects, with infinite internal motion in continuum limit, 
when the support of eigenvalues extends to the whole real axis of A. This infinite internal 
motion is the only hope to recover the perturbative QCD with the space independent 
master field. As it was mentioned in |IJ, the spectral integrals over A in the strong 
coupling expansion reproduce the lattice Feynman graphs of the large iV theory, much in 
the same way, as it took place in reduced Eguchi-Kawai models [H . The correspondence 
between the master field of our model and quenched and twisted reduced models was 
discussed in the recent paper [O]. 



1.2 Zn problem and loop averages 

These were nice features. The problems arise when one tries to introduce quarks. Within 
the large N expansion the quarks act as infinitesimal perturbation, which do not change 
the vacuum of the theory. All the physics of the conventional expansion is described 
by the Wilson loops W(C) in fundamental representation, averaged over QCD vacuum. 
The quark confinement corresponds to the global symmetry of this vacuum. The 
dynamical realization of this symmetry in QCD is the (minimal) area law W(C — > oo) ~ 
exp(— aA(C)) at large distances, combined with asymptotic freedom W(C —>■ 0) ~ 1 at 
small ones. 

It is not clear, how to get this behavior of the Wilson loop in our model. As it was 
noted already in the first paper , and discussed at length in subsequent papers || |9|, [K| , 
the effective action for the gauge field after elimination of the heavy scalars involves only 
absolute values of the Wilson loops. This leads to the spurious local Z N symmetry^] in this 
model, as well as in any other model, built from the traces in the adjoint representation 
of the gauge group. As a consequence, we get the local confinement: the Wilson loop 
vanishes, unless the loop folds on itself, so that the minimal area equals zero: 

W(C) = 5 , A (c) (1.5) 

3 Actually, we get the local U(l) symmetry, but we could always lower it to Zn by fixing the determi- 
nant of f2. 
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Formally, this corresponds to infinite string tension a = oo. 

Recently, the loop equations and the Eguchi-Kawai reduction were studied for our 
model in For the quadratic potential £/(<&), the solution of the loop equations was 
found, in which the adjoint Wilson loop, as well as the fundamental one, obey the trivial 
"zero area law" ( |1.5| ), up to O f^J corrections. The path amplitude of scalar field, 

(*|G(L)|*> = U*o ff*Wi ® ^ w *«) = A L ( tr $ 2 ) , (1.6) 

\ fc=0 / 

where the matrix <3> is treated as a vector in adjoint representation. It is implied, that 
the path {x , x±, . . . , xl} has no backtracking steps, and should it have, those steps would 
not count, as the matrices would cancel. So, this amplitude depends upon the algebraic 
length Hl4|| , L(T) of the path T, (as opposed to the usual lattice length |T| = # of steps). 
The explicit formula for ( tr $ 2 ), (in our normalization of m ), 

tr $ 2 ) = 2D ~ l (1.7) 

1 ro§(D - 1) + DJmg + 4(1 - 2£>) 



coincides with the semi-circle solution, found previously [|5] from the master field equation. 
The solution for A reads 



A = ] j\ml + (l-2D) + ^m 2 (1.8) 

These solutions, however, fail to produce any critical phenomena, which are necessary 
to match perturbative QCD. Apparently, it is impossible to induce QCD without self- 
interaction of the scalar field (if it is possible at all!). 

The usual propagator ($\G(x, y)\$) could be obtained from this one by summing over 
all paths 

(9\G{z,y)\*) = E^o 2 '^ 1 ($\G(L(T xy )) \9) (1.9) 

This sum is calculable, but the results were not presented. Presumably, they would agree 
with the solution of the corresponding wave equation || P] , inverting the Gaussian part 
of effective Lagrangean for 8p x (X). 

The same quantities could be found for arbitrary potential from the master field equa- 
tion, by noting, that the one-link integrals factorize over the path, so that the general 
invariant propagator 

(A,™ . (^nV ( ,»OL w ^-) (1.10) 
/ 1 L_1 1 
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where () n denotes the one-link averages with the measure ( |1.2|) , at fixed scalar field, 
and () $ denotes the scalar field average. As before, $ is treated as a vector in adjoint 
representation. We could write the following identity, 

^ ® n l)n j^~ y = 1 jz^ 1 M = G a(*«) (1-11) 

where V$ = ^ 



Ga(^) = (A-V*,- j P($ !B ))- 1 , (1.12) 



and 



i a In 



* * * 2D 



ip/^. (1.13) 



The details could be found in previous papers, or below, where we repeat the same 
computations for the mixed model. The function G\(Q) was computed at N = oo in || 

GM = x^m^ p (I *w -%-«» arctan ^w>) • (L14) 

with 

, x U'(v) D — l f , p(ii) 
Let us represent this function in terms of dispersion integral 

G xW = p J d „« (1.16) 
with the kernel (p\A\X) given by inverse dispersion relation 

( „ia|a> = P / (L17) 

then we see, that the one-link average acts as the operator A 



which yields the following result for the propagator 



T (A|A L |//\ (1.19) 



A 7 7 *>_0)(A-0) 

In general, there is always the trivial eigenvalue Ao = 1, corresponding to the unit matrix 
as an eigenfunction. The other eigenvalues must be less than 1, otherwise there would be 
a phase transition. 

In particular, according to H 

n xy ®nt y ) n $ y = F($ x ) (i.2o) 
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For the quadratic potential, according to ||, F(§) = A$, so that there is a trivial linear 
eigenfunction $ with the eigenvalue A. This is in complete agreement with the loop 
equations. 

The same arguments lead to the following result for the adjoint loop average 



W a (C) = 5 0M o + ^ [-1 + (l - 5 QM o) tr A L ^ 
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The leading term at N — oo, as before, arises for the backtracking loop, and it corresponds 
to the infinite string tension, regardless the solution for the scalar field density. The 
correction describes the perimeter law for the pointlike "mesons", propagating along the 
adjoint (double) path. The wave equation, effectively summing over all these paths, was 
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As for the fundamental Wilson loop, in the large N limit, we could take the SO(2N+l) 
gauge group instead of the SU (N). The large A" saddle point equations would be the same, 
up to O (jA corrections. It is also known, that the loop equations |14| of the complete 
gauge theory, with the Yang-Mills or Wilson terms, coincide with the same accuracy for 
the SU(N) and SO{2N + 1) gauge groups. 

On the other hand, there is no center in the 5*0(2^ + 1) group, hence the Wilson loop 
would not vanish. In virtue of the factorization property |14| (tr A tr B) = (tr A) (tr B) + 
0(1) the fundamental Wilson loop would be equal to the square root of the adjoint Wilson 
loop, in agreement with suggestion of 0. At A" = oo, therefore, both satisfy the zero 
area law, unless we invoke the subtle mechanism of the spontaneous generation of the 
fundamental Wilson loops in effective gauge action |7|, || |13| , which is not clear how to 
do. 



1.3 Mixed model 

At this point it is worth recalling, that our model is somewhat artificial. There was 
no physical reason in the choice of the adjoint representation of the matter field |J to 
induce QCD. We did so, simply because at that time this seemed to be the only model 
with correct counting of powers of N, which could be solved by orthogonal polynomial 
technique. Later, when the model was actually solved, the more powerful technique of 
the Rieman-Hilbert equations was found. Now, when we are no longer limited to the 
Itzykson-Zuber integral, the time came to extend the model. 

We add to the model the heavy fermion constituent field, ty x , with Nf <C A?" — > oo 
flavor components, and solve this mixed model exactly. The physical quarks q x could 
be introduced later, at larger spatial scales. The Zn symmetry, as we show, can be 
spontaneously broken in this solution, therefore the physical, light quarks would be able 
to propagate in this vacuum. 

By adjusting the parameters, we could interpolate between the local confinement and 
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free quark regimes, which gives us more chances to induce QCD. At least, the most obvious 
objection is now eliminated. 

The functional integral of the mixed model reads 



J| j D$ x exp -N tr [U{$ x )] f DV x expNti [iW^* 
n / DQ xy exp N tr «I> , < <I> V < > , , y + Vt xy ^ y H xy ^ x + O , ,/!'., //,., 



(1.22) 



<xy> 



where ^(^as) is N ® N'(N' ® N) matrix. The second dimension N' = N s Nf where Nf is 
the number of flavors, and N s = 2 Ll^J is the number of spin components. The matrices 
Hx,x+n — \ (1 + jfj) act on the spin components. 

The physical quark confinement in this case would correspond to the situation, where 
the masses of the ^fq, mesons would stay in the lattice cutoff range in the local limit, 
when the masses of physical qq mesons go to zero in lattice units. In terms of the Wilson 
loops, there would always be the perimeter law, with large decrement, due to ^ following 
the quark along the loop, but hopefully there would also be the term with area law, due 
to the induced gluon planar graphs filling the loop of the light quark. 



2 Collective fields and classical dynamics at large N 



2.1 Two densities 

The solution of the mixed model starts with the observation, that one-link integrals 



(2.1) 



depend, in virtue of the gauge invariance, only upon the densities pi, p 2 , &i, u 2 , where 

1 



P»(A) 
(ia\a x (X)\j(3) 



N 
1 

N 



tr 5(A-$,), 
tr 5(X - ® x )¥> a ¥f. 



(2.2) 



(2.3) 



with the spin indexes a, (3 and flavor indexes i,j fixed. We shall treat a as N'(&N' matrix. 

This allows one to completely eliminate \l/ from the problem, by the local change of 
variables^. Let us introduce the N'®N' matrix Lagrange multiplier e(A) for the constraint 
( [Of ), then we have to compute the integral 



J J Da J Deexp tr' (tfe($)tt) - N J dX tr'e(X)a(X) 



(2.4) 



This simple observation, applied to the old induced model, with quarks in place of immediately 
rules out the hypothesis, that it could induce QCD without the spontaneous breaking of Zn symmetry. 
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where tr' corresponds to the spin and flavor trace, and the e integration goes along 
imaginary axis. The notation e($) is used for the matrix-valued function, in practice this 
is used in the basis where $ is diagonal, where it means the diagonal matrix of e(Aj), 
where Aj are the eigenvalues of 

Integrating over \l/ first, we find 



N 



Dty exp tr' (^e($)^J oc J| det e(Aj) = exp 



i=i 



N / dXp(X) tr'lne(A) 



This yield the extra local term in effective action for the e, a variables, 
5S eff (a, e,p) = N J dX tr' (e(A)a(A) - p(A) lne(A)). 
The local field e(A) can be eliminated from equations of motion 

e(A) = p(X) HA))" 1 , 

which yields 

5S eff (a,p) = const +N J dXp(X) (-N' lnp(A) + tr'lna(A)). 



(2.5) 



(2.6) 



(2.7) 



(2.t 



The extra term in effective action for the p field was computed in ||. It starts from 
O(N) terms in the large N limit, which can be neglected in the leading order under 
consideration, provided lcJVj. Therefore, the above term in effective action serves both 
densities. 



2.2 Classical equations 

We assume, that 1 <C Nf <ti N, so that the flavor corrections are important, but still the 
classical equations for the N' <g) N' matrix a could be applied. Our objective now is to 
derive the set of equations for the one-link integral, and solve them together with these 
classical equations. For the space independent master fields p, a the classical equations 

^ 2D d 6lnI + w'm- r mu 1 ( n' p ' W + tr> a ' W ) (2 0) 

WTX6piX) +mViX) - UiX) + N{- N rtXj + tr ^ATj' (2 ' 9) 

f, 5 ^— M+^(A)MA))-\ (2.10) 

The first classical equation follows from variation of the total action with respect to 
"0(A), where Sp = tp'(X), the extra derivative being introduced to preserve the normal- 
ization condition / dXp(X) = 1. The boundary conditions are i/j(±oo) = 0. This is the 
same equation as the old one, except for the last term, coming from effective action. The 
potential 

V'(z)=[du-M- (2.11) 
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and the values at the real axis are understood as limits from the upper half plane 

V'(X + lO) = -tnp{\) + p ( dup^- (2.12) 

J A — v 

The second classical equation is new. It follows from the variation of the total action 
in o~(A), with spinor and flavor indexes implied everywhere. The mass term and the scalar 
density are proportional to the unit matrix, and the second term in general case involves 
the inverse of a matrix. We assume, that the vacuum densities p x {X) and a x {\) are 
spatially homogeneous and the second one is proportional to the unit matrix, in virtue of 
the symmetry of the model. 



2.3 Schwinger-Dyson identities 

Apart from these classical equations, which are valid only at iV — > oo, there are some 
identities, which are satisfied by the one-link integral. The simplest is the one, used in |2j 
to solve the scalar model, 

j tr (A - V^)" 1 / = tr (A - ^y 1 , (2.13) 

where 



One may rewrite this identity as follows 

tr (A - V $1 - F($ 1 ))- 1 = tr (A - Sj,)" 1 , (2.15) 

where 

F($ x ) = V^lnJ, (2.16) 
is the matrix- valued function. In the same way, as in [Q, 0L we find 

- - <™ - ™ + M- v S + tr '^) (2 - 17) 

In order to compute another variation of /, we note, that in virtue of unitarity of Q, 

dl - dl 
tr (A - V^r 1 = tr (A - ^y 1 * 2 — , (2.18) 

with sum over spins and flavors implied. On the other hand, from definition of cr(A) (for 

— = — \& ■ — — = \& (2 19) 

9* AT <5<x($)' N5a(<5>) ' 1 J 
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which yields, after moving I to the left, shifting V — » V + F, replacing the trace on the 
right by the spectral integral, 



A — /i 



with 



/ \ 1 / / \ 5 In/ 



(2.20) 



(2.21) 



From now on we are going to deal with classical solution cr(A) which has no indices in 



spin space. Using the classical equation (|2.10| ) 

N' 



2D 



(2.22) 



We cannot directly represent the left of ( |2.20[ ) in terms of the spectral integral, as 
there are derivatives V$, which prevent us from diagonalizing $. The problem is of the 
same kind, as in the pure induced model, which we solved before. We have an identical 
problem in the present model, in the scalar sector, described by (|2.15|) , so let us solve it 
first. 



2.4 First master field equation 

Let us introduce the matrix- valued function Cr A (<&) as solution of the differential equation 

V„G° A ($) = -1 + (A - F($)) G° A ($), (2.23) 
then we have the following normalization condition from Q2.15 ) 

tr ) + — ^) = (2.24) 

which holds identically for any A. We choose A to belong to the support of the eigenvalues, 
and take the real solution for G°, corresponding to the principle value prescription. 

In the large N limit, the derivative V$ acts as the following integral || 

V,G° A (*) = / dppM ^M-fW (2.25) 
J fi — <P 



which reduces the problem to the Riemann-Hilbert integral equation 

fl — V 



l + pj d^pin)^^ = (A - R(u))G° x (u), (2.26) 



m - m + - m + ^«v>, + ^ + m . ( , 27 , 



9 



For the analytic function 

T x °(z) = l + I dpp{p)^-, (2.28) 
J p — z 

which is defined in the upper halh plane and extended to the whole plane by the symmetry 
relation 

T x (z)=f x (z), (2.29) 

this is the boundary problem 

?~aV + *0) = \-R(v)+mp(v) , , 

T x °(u - i0) A - R{y) - inp(v) { ' ' 

with the well known solution 

^W= eXP (/^ arCtan A^)' < 2 ' 31 > 

The solution for G\ can be obtained from the real part of the complex conjugate boundary 
values 

7 x \v ± *0) = (A - R(y) ± wrp(i/)) = ^ ^ . (2.32) 

As long as we are interested only in the equation for density, we do not need this 
expression, but rather can use the asymptotic formula 

Jim z (l - T x °(z)) = J dpp(p)G° x (p) = MV'(X), (2.33) 



where the last relation followed from ( 2.24 ). On the other hand, we can easily find the 
same quantity from the explicit solution for T x (z) which yields the first master field 
equation 

p/ ^ + arctan _^L)= , ,,34) 
which differs from the old one by the last term in definition Q2.27 ) of R(X). 



2.5 The second master field equation 

Let us now derive the second master field equation. Repeating the same steps, we find 
for the new function 

G A ($) = (A — V$ — F($))- 1 7y($) (2.35) 
the following differential equation 

V $ G A ($) = -r ? ($) + (A - F($)) G A ($), (2.36) 

and the normalization condition 



J dpp(p)G x (p) = pj dp^f^- (2.37) 
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Furthermore, at JV = oo we arrive at the integral equation 

V (u) + pf dpp{p)^- = (A - R{y))G x {y), (2.38) 

with the same R{y) as before. The first term adds a new element to the problem. Now, 
the analytic function 



J it z — p J p — z 



in the upper half plane of z, and it is continued to the lower half plane by the symmetry, 
as before. The boundary problem 



T x {v ± iO) = (A - R{v) ± mp{v)) G x {v) ±ip [ ^^L, (2.40) 

J it v — p 

is now an inhomogeneous one. It can be solved by substitution 

T x (z) = T x °(z)J x (z), (2.41) 
which provides the following boundary problem for J x , 

(A - R{y) ± mp{v)) (j x (u ± iO)G° x (v) - G x (u)) = ±ip I (2.42) 

v ' J 7T V — p 

This equation defines the imaginary part 

G^(z/) (A — R\y)) + it l p A [y) J it v — p 
so that we could restore J x from dispersion relation 

Mz)= f ^ J ^ + '°\ (2.44) 

J It v — z 

There are no subtraction terms here, as this function should decrease as z^ 1 to satisfy 
the normalization condition. The coefficient in front of z^ 1 can be found from above 
equations; this yields the following relation 

/ d ^ Mv + ,<» = _ / %<„) + p / iv ^L (2 .45) 
The real part of this relation provides us with the second master field equation 



/ , ( pWHy) 1 (A-flQ/)) f dp r](p) ] 

p / dv \ — -q— — — — p / — = 2.46 

A — v G x (u) (A — R\v)) + it l p l [y) J it z u — p ' 



The imaginary part yields 

= / dvi){y) (2.47) 
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or, substituting the explicit formulas, 

which is quite a surprise. The "heavy quark condensate" is trivialy related to the bare 
mass. This is, actually, the lattice artifact, nothing to do with critical phenomena. It is 
convenient to renormalize a(X) — > jJm&(X), v(ty — * 2M»?(A), so that 

77(A) = p(X) - (j(A); j dvp{v) = J dva{v) = 1. (2.49) 

The expression for R(X) would not change, as it involves only the logarithmic derivative 
of <t(A), and the above equation for 77(A) would not change, being linear. 

Now, it is easy to see, that there is always the Zn symmetric solution 

a(X) = p(X) (2.50) 

of the master field equations, with the same scalar density p(X) as before. The extra term 
in -R(A) exactly vanishes in this case. 

This is the old vacuum. One may readily check, that the average = in this case. 
In general case, this average is linearly related to 77(A). As for the adjoint averages, those 
would be the same, as before, as it follows from above equations. The jf corrections drop 
from R(X), and p(A) is the same, so the heavy fermions decouple at infinite N and small 

N ■ 

However, with proper adjustment of parameters of the scalar potential, the linear 
master field equation for 77(A) (with a = p in R ) could have a nonzero solution. This is 
the spontaneous breaking of the Zn symmetry in our model. After this bifurcation point, 
the new vacuum, with two different densities would be stable. Presumably, this one does 
induce QCD, but that remains to be seen. 

Let us stress once again, that this solution does not apply to the Veneziano limit 
Nf ~ N, as in this case the classical equations for the matrix field a are no longer valid. 
However, the first correction in which we found, already breaks the Z^ symmetry of 
the vacuum. 
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4 Note added 



When this paper was finished, a new paper [15] appeared, where the trivial D = 1 
and Gaussian solutions were discussed at great length. Also, the completeness of the 
Schwinger-Dyson equations, leading to the master field equation, was questioned. Unfor- 
tunately, this objection is based on a misunderstanding. The authors of forgot about 
the gauge invariance, in virtue of which the Itzykson-Zuber integral I[§ x , &y] depend on 
tr tr but cannot depend upon tr as they suggest. The subgroup Pn of 

the gauge group independently permutes the eigenvalues <&W an d which eliminates 
the terms like J2j 

So, the gauge invariance plays the role of the "mixed" Schwinger-Dyson equations, 
they were worried about. With the gauge invariant Anzatz, the Schwinger-Dyson equa- 
tions unambiguously determine the Taylor expansion of the logarithmic derivative F of 
the Itzykson-Zuber integral, as it was discussed in ||. Actually, the uniqueness of the 
reconstruction of F is irrelevant, since we do not solve the Schwinger-Dyson identities 
for F, but rather substitute the classical equation for F into these equations, to obtain 
the master field equation for p. All we need here, is the correct Schwinger-Dyson equa- 
tion, the completeness is not used. As for the other "lessons" from the Gaussian and 
one-dimensional models, I doubt their relevance to the problem of induction of QCD. 
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